Astrophysics of magnetically collimated jets generated from laser-produced plasmas. 
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The generation of astrophysically relevant jets, from magnetically collimated, laser-produced plas- 
mas, is investigated through three-dimensional, magneto-hydrodynamic simulations. We show that 
for laser intensities I ~ 10 12 — 10 14 W cm -2 , a magnetic field in excess of ~ 0.1 MG, can collimate 
the plasma plume into a prolate cavity bounded by a shock envelope with a standing conical shock 
at its tip, which re-collimates the flow into a super magneto-sonic jet beam. This mechanism is 
equivalent to astrophysical models of hydrodynamic inertial collimation, where an isotropic wind is 
focused into a jet by a confining circumstellar torus-like envelope. The results suggest an alternative 
mechanism for a large-scale magnetic field to produce jets from wide-angle winds. 



The ejection of mass in the form of bipolar jets is ubiq- 
uitous in astrophysics, and it is widely recognized to be 
the result of a large-scale, magnetic field extracting en- 
ergy from an accreting system [T]. Models of magneto- 
hydrodynamic (MHD) jet launching P] rely on differen- 
tial rotation to shear the poloidal magnetic field compo- 
nent, B po ; = B r r + B z z, and generate a toroidal com- 
ponent, Bg, which is necessary to power and collimate 
the flow. The magnetic structure then consists of heli- 
cal field lines, with the collimation of the outflow deter- 
mined by the component of the Lorentz force perpendic- 
ular to B poh namely Fj_ = - J^V_l (rBg) + jgB po i. The 
first term of F± is the essence of self-collimation, and 
its role in jet collimation and stability has been studied 
not only through multi-dimensional simulations [3], but 
also in experiments using dense, magnetized plasmasjU- 
E]. These experiments can in fact produce flows that 
are well approximated by the Euler MHD equations, and 
whose invariant properties allow meaningful scaling of 
laboratory to astrophysical fluid dynamics [7j. However, 
collimation solely by the poloidal magnetic field, the term 
F± ~ jeBpoi, still remains to be clarified. For static 
plasma columns, the confinement was studied in the con- 
text of linear theta pinches [8]. In astrophysics, poloidal 
collimation can act in the magnetosphere-disc region on 
scales of a few au (1 astronomical unit ~ 1.5x 10 13 cm), 
where the collimation of a stellar wind depends critically 
on the the magnetic field anchored in the disc[5]. On 
scales of tens of au, it leads to the formation of axially- 
elongated cavities [TO], and may also serve to re-collimate 
potentially unstable MHD jets [IT]. On even larger scales, 
it is an essential component of the collimation of outflows 
embedded in magnetized envelopes|12j. 

In this Letter we establish the astrophysical rele- 
vance of coupling laser-produced plasmas with a strong 
magnetic field, as a platform to study jet collimation. 
Although interest in non-relativistic jet-like flows has 
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FIG. 1. Initial collimation radius (left) and time-scale (right) 
calculated with / = 0.1. Because of the dependence adopted 
for Vexp oc J 1 / 3 , the collimation time-scale is independent of 
the laser intensity. 

instigated a number of experiments using high-power 
lasers [13]. these remain so far limited to unmagnetized 
jets[]31[TS]. For laser intensities in the range 10 12 — 10 14 
W cm -2 , corresponding to laser energies E L ~ 5 — 500 J, 
with nominal laser parameters for the pulse duration of 
tt, = 1 ns, focal spot diameter of cf> = 750 /im, and wave- 
length of A = 1.064 fim, we show that under conditions 
now accessible to current facilities [16], a long-duration 
(t ^> 10 ns), strong magnetic field (> 0.1 MG) can mag- 
netically collimate jet-like flows. The basic configuration 
studied consists of a solid planar target immersed in an 
externally applied, homogeneous magnetic field Bo par- 
allel to the z-axis, and perpendicular to the target. Using 
a combination of two- and three-dimensional (3D) simu- 
lations, we provide a theoretical description of the mech- 
anism responsible for generating hydrodynamic jets via 
a conical shock, from an uncollimated plasma. These re- 
sults suggest a novel mechanism where wide-angle winds 
from stars and discs, may be re-collimated into hydrody- 
namic jets by a large-scale, poloidal magnetic field. 
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FIG. 2. (Colour online) Colour maps correspond to the logarithmic density in g cm -3 . Panels (a), (c) and (e) show a cut 
through the middle of the computational domain in the xz-plane. Contour lines in panel (a) correspond to M ma = 1 (dashed) 
and M ms = 10 (solid). Velocity vectors are shown in panel (c), while in panel (e) the contours are for the magnetic field 
intensity in MG. Panel (b) is a zoom over the conical shock region depicted in (c), and shows additionally the region where the 
flow is sub-magnetosonic, M ma < 1 (dashed line). Panel (d) is a cut perpendicular to the jet at z = 17 mm. Panel (f) shows 
the profiles on axis of density, p x 10 6 (g cm -3 ), axial velocity, v z (km/s), and ion and electron temperatures (eV). 



For a given applied magnetic field and laser energy, a 
characteristic collimation radius can be estimated from 
the equilibrium between ram and magnetic pressures 
(pv 2 ~ B 2 /8n) as E coH (cm)~ 0.8(^(J)/B (T) 2 ) 1 / 3 , 
where the bulk kinetic energy is parametrized as a frac- 
tion / of the laser energy deposited on target, Ek = fEi,. 
Numerically we find / ~ 0.3 — 0.5, which is consistent 
with experimental measurement 18J, however, consider- 
ing only the radial expansion, better estimates are ob- 
tained for/ ~ 0.1 — 0.2. The related collimation time- 
scale is estimated as t co u ~ R C oii/v e xp, where the expan- 
sion velocityQI] is v exp {cm/s)= 4.6 x 10 7 / 1 / 3 A 2 / 3 (I is 
the intensity in units of 10 14 W cm -2 and A is the laser 
wavelength in /im) . Figure [I] represents R co ii and t co u 
as a function of applied magnetic field, and shows that 
to magnetically collimate a jet-like flow with a radius of 
a few millimeters, requires field intensities > 0.1 MG, 
applied for several tens of nanoseconds. 

Numerically, we investigate the interaction of a laser- 
generated plasma plume from solid foil targets (C, Al, 
Cu) with an externally applied, steady-state magnetic 
field in the range B = 0-0.4 MG. Although in this 



regime a strong (MG) magnetic field can be generated 
from non-parallel gradients of electronic temperature 
and density [20]. it remains localized both in time and 
space[2"Tl 122). and does not affect the plasma dynamics 
over the time-scales (^> t l ) and length-scales (^> <j)) of 
interest to our work. The initial plasma evolution is mod- 
elled in axisymmetric, cylindrical geometry with the two- 
dimensional, three-temperatures, Lagrangian, radiation 
hydrodynamic code DUEDpj. coupled with SESAME 
EOS tables[24 . The plasma profiles of density, momen- 
tum and temperature (electronic and ionic) are then (at 
t = 1.2 ns) linearly mapped onto a 3D Cartesian grid 
with a superimposed uniform magnetic field, and used 
as initial condition for the 3D Eulerian, resistive MHD 
code GORGON [53 US]. We shall see that 3D calculations 
are necessary to capture the non-axisymmetric modes of 
MHD instabilities developing in the flow at late times 
(i 3> tl). Simulations were run at different resolutions 
(Ax = 35 — 65) fim and also with the initial velocity field 
randomly perturbed ( 5v/v ~ 0.05 — 0.15). The results 
are quantitatively similar, with only small differences in 
the azimuthal structure of the flow. 
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The magnetic collimation of a laser-generated plasma 
plume may be characterized by three main phases. These 
are shown in Fig. [2j for a simulation of an Aluminium 
target , with I = 1.5 x 10 14 W cm" 2 and B = 0.2 
MG. The laser propagation is anti-parallel to the z- 
axis, and the target is at z = mm. The first phase 
(Fig. [2^) corresponds to the initial expansion of the 
plasma plume, its deceleration by the radial component 
of the Lorentz force F r = jeB z , and the formation of 
a shell of shocked plasma delineating the boundaries of 
the plume. The time shown (t = 5 ns) corresponds ap- 
proximately to the maximum radial extent (R co ii ~ 3 — 4 
mm) reached by the thermally-driven expanding plasma. 
Times are given from the end of the laser pulse, unless 
otherwise stated. Because of the relatively high temper- 
atures, T e ~ 300 — 500 eV, the electrical conductivity 
is high and dissipation of magnetic flux via diffusion is 
small. This is characterized by a relatively high mag- 
netic Reynolds number, ReM = 1.4 x 10~ 20 vL/ri ~ 100, 
where v ~ 10 7 cm/s, L ~ 0.1 cm and 77 ~ 1.5 x 10~ 16 s 
are the characteristic velocity, length-scale and resistivity 
respectively. Therefore the magnetic field is "frozen" in 
the plasma, and the field lines are swept laterally by the 
flow and accumulated in the shock envelope. In addition, 
the field lines arc bent, generating a radial component 
of the magnetic field which produces an additional ax- 
ial force (F z = jgB r ). Although this is generally small 
compared to the thermal pressure gradients, we shall see 
later that the curvature of the field lines plays an im- 
portant role on the stability of the flow. Velocities in 
the plume (few xlOO km/s) are well in excess of the 
fast magneto- acoustic speed c ma , and the deceleration of 
the plasma produces a fast shock; the magneto-acoustic 
speed, c ma — \f c\ + c s) i s a combination of the Alfven 
speed, ca = B/^/Anp, and the adiabatic sound speed, 
c s = \pjpjpi with ratio of specific heats 7 = 5/3. 
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FIG. 3. (Colour online) Volume rendering of density at 25 ns, 
showing the structuring of the flow by the RT filamentation 
instability. 

The second phase is the formation of a jet via a stand- 



ing conical shock. The propagation in the axial direction 
is essentially unimpeded by the magnetic field, and the 
cavity becomes more elongated in time. The shock en- 
velope is oblique with respect to the flow velocity, and 
compresses the component of the magnetic field tangen- 
tial to the shock, while maintaining the tangential veloc- 
ity continuous. This axial focusing mechanism is eluci- 
dated in Fig. [2]:, where the velocity vectors show the flow 
being refracted across the shock, sliding along the walls 
of the cavity, before finally converging towards the axis. 
Furthermore, as demonstrated by the contour lines of the 
fast magneto-acoustic Mach number, M ma = v/c ma (Fig. 
[2^,) , the plasma in the shock envelope remains super-fast- 
magnetosonic. Its collision on axis can then generate ei- 
ther a conical shock, if the reflection is regular, or a Mach 
reflection (Fig. |2Jd) consisting of an axisymmetric triple 
shock structure, with two oblique (conical) shocks, and 
a planar shock (Mach disk) 27J. In either case, further 
acceleration of the flow and its collimation into a narrow 
jet occur as the plasma reaching the tip of the cavity is 
redirect axially by a conical shock. The whole flow con- 
figuration described so far shares many important fea- 
tures with astrophysical models of shock focused inertial 
confinement [28] . In those models the hydrodynamic col- 
limation of a (magnetically- or thermally-driven) wind 
is the consequence of the inertia of a dense, torus-like 
circumstellar envelope, which focuses the flow in the po- 
lar direction, forming prolate, wind-blown cavities, and 
jets[29- 32 . Our results show for the first time that an ax- 
ial magnetic field can in fact mimic the action of a struc- 
tured, dense envelope, and that the complex physics of 
jet collimation can be directly accessed in the laboratory. 

Finally, the third phase corresponds to the propagation 
of the jet, which undergoes one or more expansions and 
compressions that may also lead to the further genera- 
tion of shocks (interesting similarities exist with jets in 
ultra- fast accelerative flames in obstructed channels 17). 
An example of such re-focusing event can be seen in Fig. 
[2]3, where the contour lines tracing the magnitude of the 
magnetic field show a new region of compression at the 
tip of the jet (z ~ 23 mm). Figure^, which illustrates 
the plasma properties in the jet, shows the profiles on 
axis, at 26 ns, of the axial velocity, electron and ion 
temperatures, and mass density. The shock-heated jet 
has relatively low densities, and thermal equilibration 
between the ions and electrons is slow, leading to de- 
coupled temperature. The jet emerging from the conical 
shock is aligned with the magnetic field and it is poten- 
tially susceptible to the firehose instability, which may 
disrupt the flow through long (axial) wavelength, helical- 
like distortions (e.g. [33]). The condition of growth re- 
quires anisotropic pressures, P11 — P± > B 2 /An, where 
the parallel Py and perpendicular P± pressures generally 
include both the thermal pressure, and the ram pres- 
sure due to the bulk motion of the flow (pv 2 ). For the 
highly supersonic, field-aligned flows of interest here, the 
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FIG. 4. (Colour online) Line-of-sight integrated density, 
(/ pdy) in g cm -2 on a logarithmic scale, at t — 30 ns, for an 
Al target and laser intensity I = 1.5 x 10 13 W cm -2 . 

parallel pressure is Pn ~ pv 2 , and the stability condi- 
tion, assuming an isotropic thermal pressure, reduces 
to M\ — 0/3 > 1. Although this is only marginally 
met in the jet's core, the presence of a dense, strongly 
magnetized plasma at larger radii, may provide the ap- 
parent stabilization of the flow [33]. Figure [3] shows a 
three-dimensional view of the flow at 25 ns. The axial 
structure consists of alternating regions where the ra- 
dius of the flow, rf(z), and curvature of the magnetic 
field lines change from convex to concave. In the regions 
where the plasma is radially bulging out, a Rayleigh- 
Taylor type filamentation instability can develop, with 
the conditions for its growth being similar to those of a 
theta pinch •''> 1 . 35 . In particular, the growth rate, T, 
for large azimuthal mode numbers m, with wavenumber 
kg = m/rf, is given by the classical result T ~ \Jgkg, 
where g is the effective gravity at the plasma- vacuum in- 
terface, which can be approximated as g ~ v\/R c , where 
R c is the radius of curvature [35]. In addition, as the flow 
streams along the curved walls of the cavity with a ve- 
locity v > x>a, it experiences an additional centrifugal 
acceleration of the same order of magnitude. Making 
the simplifications R c ~ R co ii and v ~ v exp , which are 
valid at early times, shows that the characteristic growth 
time-scale of the instability is short, of the order of the 
collimation time-scale, tj ~ T co u/^/m. These estimates 
are consistent with the numerical results, which show az- 
imuthal perturbations, m ~ 8 — 16, growing within a few 
nanoseconds, and leading to the rapid filamentation of 
the outer edges of the cavity first, and of the jet beam 
later (see[2ji and Fig. [3]). The radially growing pertur- 
bations also propagate axially with the flow (c.f. Fig. 
[5J;) , and produce a relatively low density, broad halo sur- 
rounding the central core of the flow. 

So far we have discussed the case of a relatively strong 
magnetic field, one that is able to generate jets via a con- 
ical shock. For weaker fields in contrast, the flow stream- 
lines tend to become parallel to the magnetic field lines, 
and the flow remains instead collimated in a cylindrical 



cavity, without jets. The effects that changing the ap- 
plied magnetic field has on the collimation and morphol- 
ogy of the flow is elucidated in Fig. [2] which shows for 
a fixed laser intensity the line-of-sight, integrated mass 
density. Indeed the flow structure changes from a cylin- 
drical cavity (B < 0.03 MG), delineated by a denser 
shell of plasma, to a prolate cavity with a jet emerging 
from a focusing, conical shock. For the largest field val- 
ues (B > 0.2 MG), the focusing conical shock is closer 
to the target, and the results is a denser and narrower jet, 
which is relatively homogeneous. We note that by plac- 
ing a massive target in the jet propagation path would 
lead to the formation of a reverse shock in the jet, in a 
configuration ideal to study accretion shocks and mag- 
netized accretion columns in young stars [35]. Finally, we 
find that although different target materials lead to qual- 
itatively similar results, increasing the atomic number of 
the target, and thus the radiative losses from the plasma, 
tends to produce better collimated jets. This is a well 
known results from non-magnetized jet experiments [13]. 

The astrophysical relevance of laboratory flows rests 
on their dynamics being well approximated by ideal 
MHD 7], which implies the advective transport of mo- 
mentum, magnetic field, and thermal energy, to domi- 
nate over diffusive transport. In this regime the dimen- 
sionless Reynolds (Re), magnetic Reynolds (Rej^), and 
Peclet (Pe) numbers are much larger than unity. The 
simulations show that the bulk of the flow is well ap- 
proximated as an ideal mag neto-fluid (Re - 10 4 - 10 5 ; 
Re M ~ 100; Pe ~ 10 - 20). Thus we expect astro- 
physical simulations of related flows, performed under 
equivalent, scaled initial condition to produce qualita- 
tively similar results. From an astrophysical perspective, 
the results demonstrate that an axial magnetic field can 
on its own play the same role of a circumstellar enve- 
lope, and lead to flows similar to shock focusing models. 
Therefore providing an alternative route to explain the 
presence of jets when massive, collimating envelopes are 
not consistent with observations [37]. Moreover, the re- 
sults suggest a new framework that combines the mag- 
netic collimation of wide-angle flows with the generation 
of hydrodynamic jets, which do not suffer from poten- 
tially disruptive instabilities linked to the presence of a 
strong Bg. The predicted formation of a standing conical 
shock is also compelling, as it could possibly explain the 
presence of stationary emission features observed close to 
young stellar jet sources[38, 39J . Although the strength 
and topology of magnetic fields in those astrophysical 
jet sources remains a major open question|40). estimates 
of its intensity[41 (~ 10 mG) are consistent with those 
needed for collimation by a poloidal magnetic field. 
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